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Abstract 

We report an upper limit on the flux of relativistic monopoles based on the nonobservation of 
in-ice showers by the Radio Ice Cherenkov Experiment (RICE) at the South Pole. We obtain a 
95% C.L. limit of order 10 _18 (cm 2 s sr)^ 1 for intermediate-mass monopoles of 10 7 < 7 < 10 12 at 
the anticipated energy E to t = 10 16 GeV. This bound is over an order of magnitude stronger than 
all previously published experimental limits for this range of boost parameters 7 and exceeds 2 
orders of magnitude improvement over most of the range. We review the physics of radio detection, 
describe a Monte Carlo simulation including continuous and stochastic energy losses, and compare 
to previous experimental limits. 



1 



I. INTRODUCTION 



Magnetic monopoles are hypothetical particles carrying a net magnetic charge. Dirac[l( 
calculated the quantum unit of magnetic charge while showing that existence of magnetic 
charge leads to charge quantization. The relationship between the fundamental magnetic 
charge g and the quantum of electric charge e takes on the especially simple form g = j- in 
Gaussian units, where a is the fine structure constant. 

Magnetic monopoles occur in grand unified theories (GUTs)|2j. Most GUTs predict 
a symmetry-breaking phase transition at an energy-scale ~ 10 16 GeVQ. Such a phase 
transition can lead to localized topological defects in the form of magnetic monopoles. An 
order-of-magnitude estimate of one magnetic monopole per cubic Hubble radius gives a 
magnetic monopole number density at the GUT time of nMitnuT) ~ 10 82 m -3 , which would 
lead to magnetic monopoles dominating the Universe today [4|, in conflict with observation. 
Meanwhile, some inflationary models predict dilution of GUT monopole density to hopelessly 
undetectable levels^. 

No monopoles have been verifiably detected to date[H|. Reports of magnetic monopole 
detections 0, 0, @] have been challenged, often by the original authors themselves (lol. 11 . 
Alternate explanations range from ionized nuclei [uj to hardware malfunctions 0] to new 
physics other than monopoles 13 • 

In this paper we present limits on relativistic monopole fluxes determined from five years 
of data collection by the RICE experiment [12]. As noted by Wick et al. [13] the RICE array 
is particularly well-suited to ultrarelativistic monopole detection because of a combination of 
large effective volume and favorable scaling with energy. Our limits for fluxes over the range 
of monopole Lorentz boost parameters 10 8 < 7 < 10 12 are about 100 times more restrictive 
than the tightest published limits from AMANDA 0, Baikal[l5|, and MACRO [l6|, and 
more than 1000 times more restrictive than the original Parker bound. 

This paper first reviews monopole properties and previous limits in Sec. [TTJ We then 
provide an overview of RICE in Sec. IHIl Sections IIVI IVj and [VI] discuss monopole energy 
loss, our Monte Carlo simulation thereof, and RICE's response. Results are summarized in 
Sec. [Vnl 



II. LIGHT RELATIVISTIC MONOPOLES 
A. Formation 

Until monopoles are detected their masses should be considered unknown. Although 
GUT-scale monopoles are commonly believed to be extremely heavy (10 17 GeV) and unde- 
tectable rare as a result of inflation, there are other mechanisms resulting in production of 
much lighter monopoles after inflation. Kephart and Shafi[17] proposed that later symmetry- 
breaking events at lower energies could occur, resulting in monopoles with a magnetic charge 
of a few Dirac units and masses in the range 10 7 — 10 13 GeV. Other candidate GUTs lower 
the GUT energy scale, allowing for 10 8 GeV monopole mass ranges, and conformal field 



theories can be developed allowing for monopole masses on the order of 10 5 GeV 13j. Be- 
cause magnetic charge is conserved, magnetic monopoles formed in the early universe should 
persist through the current epoch[l8[]. In fact, these lighter monopoles have been suggested 
as possible ultra-high energy cosmic rays beyond the GZK[l9| cutoff [l~3| . 



Such considerations motivate the search for so-called intermediate-mass monopoles 
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(IMM's), that is, monopoles with mass significantly less than the conventional GUT en- 
ergy. Our study concentrates on this mass range because such monopoles are expected to 
be ultrarelativistic and readily detectable by the RICE array 



B. Experimental Limits To Date 

In 1970, Parker pointed out that the abundance of magnetic monopoles is constrained 
by the requirement that magnetic monopole currents be insufficient to deplete the galactic 



magnetic field [20j. Only in the past decade have astrophysical experiments been able to im- 
prove upon the original Parker bound of flux ~ 10~ 15 (cm 2 s sr) _1 [16|. The first observational 
astrophysics experiment to obtain limits stronger than the Parker bound was MACRO (the 
Monopole Astrophysics and Cosmic Ray Observatory) at Gran Sasso, Italy. The MACRO 
bound for monopole velocities v = /3c > 0.99c is a flux upper limit of 1.5 x 10 _16 (cm 2 s sr) _1 . 
Upper bounds of this order of magnitude were also reported for 4 x 10~ 5 < (3 < 0.99(l6|. 

Since MACRO'S termination in 2000, "neutrino telescopes" have conducted searches for 
relativistic magnetic monopoles. Both AMANDA (the Antarctic Muon and Neutrino Detec- 
tor) and the Baikal Neutrino Telescope have reported flux limits stronger than MACRO'S 
for (3 > 0.8. Baikal has set an upper bound on monopole flux of 5 x 10~ 17 (cm 2 s sr) _1 for 
13 ~ 1 and 10 7 < M < 10 14 [l5|, [2l|. AMANDA's most recent limit for monopole masses 
10 8 < M < 10 11 is 2.8 x 10 _17 (cm 2 s sr)" 1 at (3 « lQ- Preliminary results from IceCube's 
9-string configuration indicate that new lower bounds [somewhat below 10~ 17 (cm 2 s sr) _1 at 
(3 ~ 1] may be forthcoming [22| . 

The SLIM (Search for Light Magnetic Monopoles) experiment at the Chacaltaya High Al- 
titude Laborary in the mountains of Boliva offers complementary sensitivity to IMM's. This 
nuclear track detector experiment is designed especially to search for light monopoles (mass 
10 5 GeV to 10 12 GeV) over a wide range of velocities (including (3 > 4 x 10 -5 for 1-Dirac- 
charge monopoles). SLIM's latest monopole flux limit at (3 « 1 is 6.5 x 10~ 16 (cm 2 s sr) _1 for 
Earth-crossing monopoles, or 1.3 x 10~ 15 (cm 2 s sr) -1 if upgoing monopoles are blocked[23]. 

Stronger flux limits based on astrophysical considerations [such as an "extended Parker 
bound" of less than 3 x 10 _22 (cm 2 s sr) _1 for IMM's, based on a more realistic model of 
galactic magnetic fields [24]] have also been proposed. 



C. Relativistic IMM's 

Because of their moderate mass, IMM's may acquire highly relativistic velocities. Wick 
et al. [13] use a model of monopole traversal of intergalactic magnetic fields (similar to the 
model underlying the Parker bound) to estimate that IMM's created in the early universe 
would now have typical kinetic energies on the order of 10 16 GeV, with a comparable spread 
in energy. PeV-mass monopoles would therefore reach boost factors 7 ~ 10 10 . The fact that 
IMM's acquire such "ultrarelativistic" 7 values provides a mechanism for their detection. 
Any particle travelling through a medium loses energy, but ultrarelativistic charged particles 



do so dramatically, initiating showers in the surrounding medium[25j. It is through detection 
of such showers that RICE is sensitive to magnetic monopoles. 
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III. RICE 



RICE (the Radio Ice Cherenkov Experiment) consists of 16 data-taking antennas buried 
in the Antarctic ice at the Martin A. Pomerantz Observatory (MAPO) about 1 km from the 
geographic South Pole. The antennas, which are roughly contained within a cube of ice ~200 
m on a side with its center ~ 150 m below the surface, have peak sensitivity in the 200-500 
MHz regime. Collisions of highly energetic neutrinos (E > 10 s GeV) with in-ice atomic 
nuclei result in electromagnetic and/or hadronic cascades, depending on neutrino flavor [12]. 
The superluminal velocity of the charged particles comprising the cascade creates coherent 
broadband Cherenkov radiation in the radio frequency domain. The broad frequency spec- 
trum leads to a sharp detectable pulse of Cherenkov radiation 12|, 26|, 27] in the time domain. 



The pulse propagates outward from the shower axis in a cone with opening angle 8 C given 
by B 

* e = arccos|±|. (11 



Here n is the index of refraction of the medium evaluated at the radiation frequency. For 
highly relativistic particles (j3 ~ 1) in ice (n ~ 1-78), the Cherenkov angle is about 0.97 rad. 

A. RICE hardware and data taking 

We now briefly summarize the RICE operation, referring the reader to more detailed 
descriptions found elsewhere 12| . 



When the RICE detector is live, a trigger occurs if four or more antennas register high- 
amplitude voltages within a time coincidence of 1.25 microseconds. Triggers initiate an 8.192 
microsecond waveform capture, sampled at a rate of 1 x 10 9 samples per second, for all under- 
ice antennas. Wave forms are approximately centered on the time that the trigger latched 
the data acquisition system. Software surface vetoes (maximum rate ~80 Hz) and hardware 
surface vetoes (maximum rate ~200 kHz) suppress anthropogenic backgrounds originating 
above the array on the surface. Offline, the signal profiles are fed through a reconstruction 
algorithm, described in[28], which is designed to reject spurious triggers while preserving 



a large fraction of high-energy cascade signals. During data collection from 1999 to 2005, 
RICE was operated for a livetime of 74.1 x 10 6 s, resulting in 1.035 x 10 6 triggers. All triggered 
events were ruled out as viable neutrino candidate events during offline analysis, from which 
an upper bound on the diffuse high-energy neutrino flux was derived|28[]. Because relativistic 
monopoles also generate in-ice cascades, this same data set can be used to derive an upper 
bound on the flux of relativistic monopoles. To produce a bound it is necessary to quantify 
the sensitivity of RICE to monopoles. A Monte Carlo simulation code ("monoMC") was 
therefore created to evaluate RICE's monopole detection efficiency, as detailed below. 

IV. MONOPOLE ENERGY LOSS IN MATTER 

Our model of monopole energy loss is based on the muon/tau energy loss model of Dutta 
et al. [29]. In this model, energy loss by a muon traversing a medium is expressed as 

dE 

- = a + ( JE. (2) 
ax 
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The a term is the energy loss per distance (units: g/cm 2 ) due to ionization of the medium. 
The (3 term [30] is the sum of three terms reflecting bremsstrahlung, pair production, and 
photonuclear effect energy losses. Each energy loss mechanism is calculated separately. 
Defining y to be the fraction of its energy lost by the particle in a single interaction with 
the medium, each of the three terms $ is found by convolving y with the partial interaction 
cross section with respect to y <^fy: 

l3i{E) = -j / y dy (3) 

Here, N is Avogadro's number and A is the average atomic mass number of the medium. 
The full formulae for a and the ?/ m i ni , Z/max,, and da^/dy needed to calculate each $ are 
given in [29]. In general, the expressions are functions of particle mass and energy and of 
various properties of the medium, although a and the individual /3j's are only weakly energy 
dependent. 



A. Discrete Loss Processes 

Although energy loss due to ionization can be treated as smooth and continuous with 
little loss of accuracy, we explicitly model the stochastic fluctuation in pair production and 
photonuclear energy losses. Combining Eq. (j2J) with and replacing the integral in 
with the corresponding Riemann sum gives the result (jl]), where AEi is the energy loss via 
process % (brem., pair, or photonuclear) over a small distance Ax: 

Recasting the energy loss equation this way effectively sorts the total energy loss into an 
arbitrary number of bins, each of which spans a length Ay of the possible y values. Since 
yj is the fractional energy loss in a single interaction within bin j and E is the total energy 
of the particle, (yjE) is the energy loss for a single interaction in the j th bin. Each term 
of the Riemann sum represents an energy loss, so if (yjE) is the energy lost in a single 
interaction, the remaining multiplicative terms in the summand give the expectation number 
of interactions in the jth bin (rijj): 

N ddi . 

(n„) = -j Ax-A, (5) 

Therefore, accurately modeling the stochastic variation in bremsstrahlung, pair production, 
and the photonuclear effect is equivalent to replacing (riy) by a random number drawn from 
a Poisson distribution of expectation value (riij) when evaluating the energy loss expressions 



B. Generalization to Monopoles 

Only a few changes are needed to convert the stochastic model of muon energy loss to 
a model of magnetic monopole energy loss. First, the muon mass must be replaced by the 
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magnetic monopole mass. Because bremsstrahlung falls off by inverse powers of particle mass, 
the bremsstrahlung energy loss contribution is negligible for even light magnetic monopoles 
and will be subsequently disregarded [1 31] . It should be noted that at large masses (^1 TeV), 
/3 P air production can become difficult to calculate numerically due to rounding error. However, 
pair production energy loss approaches an asymptotic limit with increasing particle mass 
and varies with mass by only a few percent for masses above ~ 100 MeV. 

Next, due to Dirac's quantization condition a magnetic monopole of 1 Dirac charge will 



lose energy equivalent to an electric charge of l/(2a) times the proton charge [251] . Accounting 
for this large effective charge only requires multiplying the expectation number of interactions 
by l/(2a) 2 « 4700. 

This procedure for modelling magnetic monopole energy loss in matter has assumed the 
particle to be a simple Dirac monopole, that is, a point source of magnetic charge with 
no further internal structure. "Actual" magnetic monopoles may contain internal color 
fields and lose additional energy through hadronic interactions beyond the photonuclear 
effect [3ll] . However calculations of such processes are highly model-dependent (l3| and not 
further considered in this analysis. 

Figure [T] shows the average monopole energy loss in ice and "standard rock" (A=22, 
p = 2.65 g / cm 3 (29| ) . along with the same results for muons. While the muon mass is fixed, 
the monopole rest mass is constrained to vary inversely with gamma such that total energy 
is fixed at a reference energy of 10 16 GeV. Figure [2] shows various contributions to the energy 
loss of a 10 16 GeV monopole. Much of the difference between muons and monopole energy 
losses is due to the monopole's large effective charge. The curves indicate average energy 
loss due to the three principal mechanisms, while the points show actual stochastic energy 
loss (as averaged over a 50 m interval). The photonuclear effect is the dominant energy loss 



mechanism at 7 > 10 , while ionization energy losses dominate below this value [321] . Because 
the photonuclear mechanism results in hadronic showers generated by nuclear recoils, we may 
ignore LPM[33| effects. 

There are substantial uncertainties in extrapolations of photonuclear losses to ultrahigh 
energies. The primary unknown is the hadronic contribution of real photon-nucleon scatter- 



ing. Consider an extrapolation based on the Froissart bound setting in at about 50 GeV [34 
The energy dependence of the photon-nucleon cross section u^n of such a model is 

ct 7 jv(£ 7 ) = 114.3 + 1.671n 2 (0.0213£/GeV) /xb (6) 

This 1981 model of Bezrukov and Bugaev 35[ predates the discovery at HERA of parton 
distributions at small— x that causes cross sections to grow like fractional powers. For com- 
parison, the 1998-2001 post-HERA photonuclear cross section of Donnachie and Landshoff 
[361 ] is within 10% of Eqn. at = 10 6 GeV, while being about 16 times larger at 10 11 
GeV. As shown in Fig. 3 of Ref. [29], the more complicated cross sections developed by 
extrapolating structure functions track the simple expression of Eq. [6] very well. Using small 
and slow-growing cross section models such as these is in some sense conservative. It results 
in dim showers that are less likely to trigger. However, small cross sections also tend to 
develop fewer showers failing the time-over-threshold cut discussed below, and fewer showers 
are stopped by the Earth. No one cross section model is in all cases the "most conservative." 
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FIG. 1: (Color online) Muon energy loss in ice and standard rock compared to monopole energy 
loss in ice and standard rock, as a function of boost parameter 7. Monopole results are shown for 
total energy 10 16 GeV (the energy assumed in this analysis) and, for comparison, 10 14 GeV. The 
figure shows that, over the kinematic range of interest, monopole energy loss depends strongly on 
7, but for a given 7 it is only weakly mass dependent. 




FIG. 2: (Color online) Total energy loss versus 7 for 10 16 GeV monopoles, showing stochastic 
variation over 50m intervals. Lines show average contributions from different processes. Also 
shown, again for comparison, is the average total energy loss for 10 14 GeV monopoles. In the 
latter case, stochastic variation and energy loss contributions from the various processes are almost 
identical to the former case. 



V. MONOPOLE MONTE CARLO 

The monopole energy loss model described above has been incorporated into a Monte 
Carlo simulation of a magnetic monopole travelling near RICE. The Monte Carlo simulation 
randomly generates a monopole trajectory and energy loss, then determines the voltage 
response of the RICE antennas. Each voltage profile is subsequently processed using the 
same reconstruction software as was used in the 2006 RICE data analysis [2^]. 

With the Monte Carlo simulation, we generate monopoles at the reference energy (10 16 
GeV) with seven 7 values logarithmically spaced from 10 6 to 10 12 . The maximum value of 
7 = 10 12 corresponds to a 10 TeV monopole mass. Smaller masses are probably unrealistic 
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FIG. 3: (Color online) Final energy of magnetic monopoles with initial energy 10 16 GeV after 
crossing the Earth. The angle labeling each curve is the opening angle between the monopole's 
velocity vector and the zenith at the monopole's point of exit. 

theoretically. The lower 7 bound is somewhat ad hoc, chosen to bracket the kinematic regime 
to which RICE is most sensitive. 

For each simulated monopole, the Monte Carlo first generates a random monopole tra- 
jectory with an impact parameter (distance of closest approach to the array) less than a 
7-dependent maximum impact parameter r m [37|. The r m value used for a given simulation 
series must be large enough to include virtually all detectable monopoles without being so 
large that zero or a negligible number of the simulated monopoles are detected, in which 
case no flux bound can be calculated (or the flux bound will be unnecessarily high after ac- 
counting for statistical error). The selected r m values were based on preliminary simulations 
and are shown in Table [H 

A. Passage through Earth 

Unlike ultra-high energy neutrinos, not all upcoming monopoles range out before reaching 
the detector, although their energy loss in-transit can be substantial. The terrestrial density 
integrated over the length of the monopole's path |g/c m 2 |, as a function of approach angle, 
is taken from the Preliminary Reference Earth Model [38j]. The monopole energy loss in 
transiting this material is calculated over 50 increments of equal column thickness. The 
calculation assumes A = 22 for material within the Earth; however, the photonuclear effect, 
which dominates energy loss, is independent of A aside from a nuclear shadowing factor 
affecting the result at the level of a few percent. If at any time the monopole's 7 falls below 
2, the monopole is considered to have effectively stopped. In practice much larger 7 values 
are needed for RICE to trigger. 

Figure [3] shows the energy remaining after crossing the Earth for monopoles initially 
having energy 10 16 GeV with a range of incidence angles. Energy loss increases with 7, so 
monopoles with 7 ;$ 10 7 lose a negligible fraction of their incident energy. At higher 7, 
energy loss can be significant. Beyond 7 ^ 10 10 , the Earth is opaque to monopoles. 
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B. Propagation Through Ice 



After calculating the energy lost by a monopole en route to the ice target, the monopole's 
interaction with the ice itself is simulated. So long as the monopole is "out of range" of the 
RICE antennas, it is propagated along in 10 m steps, with its 7 value decremented in each 
step to account for energy lost in travelling that 10 m path. A monopole is considered to be 
"in range" as long as at least one of the RICE antennas is within 0.33 rad of the Cherenkov 
cone anywhere along a 900 m segment extending from 550 m ahead of the monopole to 350 m 
behind it. (We choose an extended path length to accurately model the time-over-threshold 
cut in RICE's offline analysis, described later.) The angle 0.33 rad is 1.70 half widths of 
the Cherenkov radiation's angular distribution at 0.2 GHz(26l. 39| . which is the lower edge 



of RICE's frequency sensitivity. At higher frequencies the Cherenkov radiation's angular 
distribution is more tightly confined. The path cutoff distances of 550 m and 350 m were 
chosen to guarantee that all signal arriving at an antenna in the first 1.2 /xs following the 
Cherenkov peak will be considered "in range" [40] . 

Once a monopole comes into range, the size of the simulation steps is reduced from 10 
m to 0.4 m. As a result of this small step size, all the energy lost within the interval can 
be treated as originating at a single "subshower" with a pointlike source, while introducing 
signal arrival timing errors no greater than 0.4 ns. By comparison, the actual experiment's 
digitizer samples every nanosecond. 

The complete voltage profile V(t) at each antenna is determined by coherently summing 
the voltage contributions of various subshowers in each time bin. Since the photonuclear 
effect dominates energy loss, the initial energy of each subshower is taken to be equal to the 
monopole's total energy loss within the corresponding distance interval. The signal phase 



will vary slowly with viewing angle [41|]; however, in most cases the dominant emission arises 
from a coherent region along the track centered around the Cherenkov point. This region 
is typically of length y\R ~ 10 — 50 m and subtends a small viewing angle. In rare cases, 
single large subshowers at modest distances dominate the radiation, but in this case the 
absolute phase of a single subshower is irrelevant. Accordingly, we ignore all signal phases 
other than from travel time in our analysis. 



VI. RICE RESPONSE 
A. Signal Analysis 

RadioMC, a Monte Carlo simulation of the radiofrequency signals caused by cascades 



(discussed in Ref. [281]) has been previously developed in the context of RICE's high-energy 
neutrino flux studies. Because of the necessity of modeling a large number of subshowers for 
each monopole, we use a streamlined version of RadioMC that does not take into account 
the depth-dependence of the index of refraction in the firn, defined as the upper 175 m of 
ice. This introduces two errors for which we correct as follows: First, in order to account for 
the presence of a "shadow zone" j28| , voltage contributions are suppressed when a monopole 
is more distant than a depth-dependent limiting horizon. Although this horizon increases 
with both source depth and antenna depth, it blocks all but two of the antennas from 
seeing even the deepest events at distances greater than 20km. Second, the streamlined code 
overestimates signal arrival time by assuming that firn ice has the same index of refraction as 
deep ice (n = 1.78), whereas n actually falls off to 1.37 at the surface. We estimate the timing 
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FIG. 4: (Color online) Path geometry of a nearby, but otherwise typical, downgoing monopole 
(7 = 10 7 ). See explanation in text. 



error from each subshower's signal using a parameterization of index of refraction vs depth, 
and time-shift each subshower's voltage contribution accordingly. This effect translates the 
Cherenkov point along the path, as a first-order correction for refractive effects. 

From the voltage profiles generated by RadioMC, we determine which monopoles would 
cause the array to trigger. RICE triggers when four or more antennas exceed a voltage 
threshold of roughly 200-400 mV (after amplification, and depending on the background 
levels at a given time) within 1.25 [is. As the precise threshold varied over the life of the 
experiment, the threshold used in any given Monte Carlo simulation is randomly selected 
from the historical threshold distribution. A typical example of a "trigger" is shown in Figs. 
H] and In Fig. HJ the long, diagonal line indicates the downgoing monopole's path, and 
the cones show the Cherenkov radiation emitted at the boundaries of the "in range" path 
segment. The missing section of the first Cherenkov cone (at upper-left) indicates where the 
cone intersects the ice surface. The square corresponds to the surface area mapped in Fig. 



3 of[12|, with the RICE antennas (enlarged for clarity) arranged in the ice below. For scale, 
the MAPO building shown on the surface is ~50 m long. Figure [5] shows the voltage signal 
generated by the same monopole. 

As a final analysis step, following the procedure of [28], each calculated monopole voltage 
profile is embedded in an unbiased event (a random recording from the actual RICE an- 
tennas, including thermal noise, etc.). The combined signal is fed through the same event 
reconstruction algorithm used in the 2006 neutrino analysis, to determine which monopoles 
would survive the offline event reconstruction analysis. 



B. Flux Upper Bounds 

Supposing that a number n from a sample of N monopoles trigger the array and survive 
reconstruction in a simulated sample with maximum impact parameter r m , an upper bound 
on monopole flux is derived as follows. Although the literal detection efficiency is e t — 
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FIG. 5: Simulated voltage vs time (as measured at the data acquisition system) in each RICE 
antenna channel as caused by the monopole of Fig. 01 The voltage graphs for the different antennas 
have been shifted vertically for clarity. 



we take the efficiency to be the lower bound of a 90% Agresti-Coull confidence interval [4^ 
about e t . Then, e s has only a 5% chance of being greater than the efficiency that would be 
obtained in an infinite number of trials. The cross section for monopole detection is 



a 



7rr„ 



i e sfb 



(7) 



Here fb is a factor which reflects the estimated decrease in detection due to the birefringence 



effect, which is not directly modeled in the Monte Carlo 43|. The 0.12% birefringence implied 



by measurements near Taylor Dome correspond to a sensitivity reduction of approximately 
14% that is roughly independent of energy; recent measurements at the South Pole indicate 
no observable birefringence (<0.01%) to a depth of ~1.1 km 44]. We conservatively take 
/fe=0.86 for the purposes of calculating our upper limit. 

Assuming Poisson statistics, the 95% upper bound on monopole flux is 

2.995 



4vrL 



a 



The factor 4tt converts the flux from units of (cm 2 s)~ x to (cm 2 s sr) _1 , assuming an isotropic 
distribution of monopoles; L is the livetime of the experiment. Although the 2006 RICE 
neutrino analysis incorporated data from 1999 through 2005, the 1999 and 2000 data are not 
considered for this monopole analysis due to differences in the detector configuration between 
those earlier data and subsequent (>2000) datasets. Thus we use a livetime of 58.3 x 10 6 s, 
corresponding to the 2001-2005 RICE operations, for calculation of our sensitivity. 



VII. RESULTS 

Figure M shows the upper bound on magnetic monopole flux as a function of 7 for 
monopoles of initial energy 10 16 GeV. These limits are based on a few to a few dozen 
simulated detections per 7 out of a generated sample of ten thousand monopoles per 7, as 
tabulated in Tabled [45j. Our limits have been degraded by the statistical uncertainty in 
our Monte Carlo-estimated efficiency for each 7 bin. 
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FIG. 6: (Color online) Upper bound on magnetic monopole flux from RICE's null observation 
in 2001-2005 versus monopole mass an d 7 . Also shown are the Parker bound and results from 
MACRO 0, AMANDaQ], and Baikal[l5]. 
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TABLE I: Final Monte Carlo results as a function of monopole 7 and mass (columns 1 and 2). 
Subsequent columns show Monte Carlo parameters (number of simulations, maximum impact pa- 
rameter r m , number of monopoles passing software reconstruction, and the resulting cross section 
for monopole detection by RICE). The resulting bound on monopole flux is tabulated in column 7 
and plotted in Fig. [6J 



We note that across a range spanning 4 orders of magnitude in 7, the flux of 10 GeV, 
single-Dirac-charge monopoles incident on the Earth is less than lxl0 _18 (cm 2 s sr) -1 . Below 
7 fa 10 8 , the utility of RICE for IMM detection rapidly deteriorates. 

For a given initial energy (as assumed here), several factors contribute to the 7 dependence 
of the upper bound on monopole flux. Monopole energy loss, and hence signal strength, 
rises with 7. Below 7 ~ 10 7 , few monopoles are bright enough to cause the array to trigger. 
Sensitivity peaks at 7 « 10 s and falls off again at higher 7. The blocking of upgoing 
monopoles for 7 > 10 10 contributes to this high-7 sensitivity falloff. However, a more 
significant effect at 7 ^ 10 9 is decreasing efficiency as more monopoles fail the "time-over- 
threshold" cut. This cut is a requirement in the offline analysis requiring that the signal 
voltage not persist for more than 800ns after initially exceeding threshold. It serves to filter 
out long-lived anthropogenic noise sources. Because of this criterion, only 49% of triggering 
monopoles are rejected in the offline analysis at 7 = 10 7 , but fully 94% are rejected at 

7 = 10 12 . 
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0.57 



TABLE II: Error analysis simulations, the number of detections in each (out of 10000 at 7 = 10 10 
using r m = 20 km), and the consequent cross section for monopole detection. 

A. Systematic Errors 

Preliminary simulations were used to select maximum impact parameters r m such that 
the effect of ignoring more distant monopoles would be negligible. Also, because the value e s 
is used in place of the raw efficiency e t , our limits are less prone to being artificially tightened 
by statistical fluctuations (although this conservative choice does weaken them by as much 
as 60% in the case of 7 = 10 12 ). 

Inaccurate modelling of voltage amplification, radio signal attenuation in ice, signal phase, 
monopole energy loss, and initial monopole energy all provide potential contributions to 
systematic error. In order to quantify these errors, the simulation was modified and rerun 
5 times at the intermediate 7 = 10 10 . A summary of our systematic checks is presented in 



First, we reanalyzed our sample of 10k monopole energy loss profiles at 7 = 10 10 , modi- 
fying each of several RadioMC parameters one by one. 

Reducing the RICE amplifiers' voltage amplification by a factor of 2 in each channel does 
not cause a statistically-significant (at the 95% level) change in detection sensitivity. This is 
expected because most detected monopoles produce signals well in excess of the triggering 
thresholds. 

Reducing the in-ice signal attenuation length by a factor of 2 at all frequencies also 
causes a statistically insignificant change in sensitivity. Increasing signal attenuation causes 
some otherwise detectable monopoles to become too faint to trigger, but it also allows some 
monopoles that would otherwise fail the time-over-threshold cut to pass. These competing 
effects largely cancel each other out. 

To explore the effects of variation in signal phase, we consider the results of subjecting each 

subshower's contribution to a phase shift of (%) ( — ) , where A6 is the difference between the 



viewing angle and Cherenkov angle and oq is the frequency-dependent half-width (divided 
by a factor of 1.17) of the Cherenkov radiation's angular distribution j39|. This is a large, 
ad-hoc phase variation that we used to test the potential effects of phase shifts on sensitivity. 
Again, no statistically significant change in detector sensitivity occurs. This indicates that 
our results are relatively insensitive to phase variations with viewing angle. 

Next, we generated a new sample of 10k monopole energy loss profiles, assuming a factor of 
10 increase in monopole energy loss. This extreme change causes a modest (and statistically 
insignificant at the 95% level) decrease in monopole detections. From Fig. [21 it can be 
seen that the results of raising energy loss by a factor of 10 will be similar to the results of 
increasing 7 by a factor of about 8, so the small magnitude of this change is expected in 



Table HD 
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light of the earlier results (Fig. [H]). 

Finally, we generated a second new sample of 10k monopole energy loss profiles, assuming 
a factor of 100 decrease in initial monopole energy. This lower initial energy is more con- 
sistent with, for example, the initial energy assumed by the Baikal Collaboration for their 
analysis jl5j. The resulting sensitivity decrease is statistically insignificant. Lowering total 
energy at fixed 7 is equivalent to decreasing the monopole's rest mass, and energy loss is 
only weakly mass dependent (Fig. [1]). While decreasing initial energy may sometimes cause 
more monopoles to be trapped in the Earth, this effect makes no difference at 7 = 10 10 , for 
which the Earth is opaque to monopoles in either case. 

In the above discussion, we have assumed a much greater uncertainty in energy loss and 
initial energy than in signal attenuation or amplification. In general, RICE's sensitivity to 
relativistic monopoles is largely determined by the detector/ice geometry, with details such 
as energy loss and signal phase shift models playing much smaller roles. The photonuclear 
component of monopole energy loss is the least well-established aspect of the energy loss 
model. Although variations by an order of magnitude are seen in the literature (Sec. HVBj) . 
these variations are unlikely in any case to alter our limits by more than a factor of 2. 



VIII. SUMMARY 

From the nonobservation of highly ionizing shower "trails" we have derived the monopole 
flux upper limits shown in Fig. [61 which are on the order of 10~ 18 (cm 2 s sr) _1 . Previously, 
AMANDA[li], Baikalfll], and MACRO [l6| determined monopole flux limits on the order 
of 10~ 16 (cm 2 s sr) _1 for (3 greater than 0.8, 0.8, and 4 x 10~ 5 , respectively. Although the 
results of this study cover a much narrower range of (3 values than previous works, it is the 
range that is of the greatest interest for IMM searches. Within much of this kinematic range 
(E = 10 16 GeV; 7 > 10 8 ), monopole flux limits from RICE are stronger than the limits from 
any previous astrophysical monopole search by more than an order of magnitude. 
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